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Abstract 

The dynamics of phase transitions plays a crucial role in the so- 
called interface between high energy particle physics and cosmol- 
ogy. Many of the interesting results generated during the last 
fifteen years or so rely on simplified assumptions concerning the 
complex mechanisms typical of nonequilibrium field theories. Af- 
ter reviewing well-known results concerning the dynamics of first 
and second order phase transitions, I argue that much is yet to 
be understood, in particular in situations where homogeneous 
nucleation theory does not apply. I present a method to deal 
with departures from homogeneous nucleation, and compare its 
efficacy with numerical simulations. Finally, I discuss the inter- 
esting problem of matching numerical simulations of stochastic 
field theories with continuum models. 
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LECTURE I 



1 Homogeneous Nucleation 

The fact that the gauge symmetries describing particle interactions can be 
restored at high enough temperatures has led, during the past 15 years or so, 
to an active research program on the possible implications that this symme- 
try restoration might have had to the physics of the very early Universe. One 
of the most interesting and popular possibilities is that during its expansion 
the Universe underwent a series of phase transitions, as some higher symme- 
try group was successively broken into products of smaller groups, up to the 
present standard model described by the product SU (3)c <8> SU (2) l <8> U(l)y. 
Most models of inflation and the formation of topological (and nontopolog- 
ical) defects are well-known consequences of taking the existence of cosmo- 
logical phase transitions seriously [[!]]. 

One, but certainly not the only, motivation of the works addressed in this 
talk comes from the possibility that the baryon asymmetry of the Universe 
could have been dynamically generated during a first order electroweak phase 
transition ||. As is by now clear, a realistic calculation of the net baryon 
number produced during the transition is a formidable challenge. We prob- 
ably must invoke physics beyond the standard model (an exciting prospect 
for most people), push perturbation theory to its limits (and beyond, due 
to the nonperturbative nature of magnetic plasma masses that regulate the 
perturbative expansion in the symmetric phase), and we must deal with 
nonequilibrium aspects of the phase transition. Here I will focus on the lat- 
ter problem, as it seems to me to be the least discussed of the pillars on 
which most baryon number calculations are built upon. To be more specific, 
it is possible to separate the nonequilibrium aspects of the phase transition 
into two main subdivisions. If the transition proceeds by bubble nucleation, 
we can study the propagation of bubbles in the hot plasma and the transport 
properties of particles through the bubble wall. A considerable amount of 
work has been devoted to this issue, and the reader can consult the works of 
Ref. for details. These works assume that homogeneous nucleation the- 
ory is adequate to investigate the evolution of the phase transition, at least 
for the range of parameters of interest in the particular model being used to 
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generate the baryon asymmetry. This brings us to the second important as- 
pect of the nonequilibrium dynamics of first order phase transitions, namely 
the validity of homogeneous nucleation theory to describe the approach to 
equilibrium. This is the issue addressed in this talk. 

Nucleation theory is a well-studied, but far from exhausted, subject. Since 
the pioneering work of Becker and Doring on the nucleation of droplets in 
supercooled vapor Q, the study of first order phase transitions has been 
of interest to investigators in several fields, from meteorology and materi- 
als science to quantum field theory and cosmology. Phenomenological field 
theories were developed by Cahn and Hilliard and by Langer || |6| in the con- 
text of coarse-grained time-dependent Ginzburg-Landau models, in which an 
expression for the decay rate per unit volume was obtained by assuming a 
steady-state probability current flowing through the saddle-point of the free- 
energy functional || |7j|. The application of metastable decay to quantum 
field theory was initiated by Voloshin, Kobzarev, and Okun ||, and soon 
after put onto firmer theoretical ground by Coleman and Callan |§. The 
generalization of these results for finite temperature field theory was first 



studied by Linde |fL0|l , and has been the focus of much recent attention |TT 

The crucial ingredient in the evaluation of the decay rate is the compu- 
tation of the imaginary part of the free energy. As shown by Langer 0, the 
decay rate 1Z is proportional to the imaginary part of the free energy JF, 

n = , (i) 

where E_ is the negative eigenvalue related to metastability, which depends 
on nonequilibrium aspects of the dynamics, such as the growth rate of the 
critical bubble. Since JF = —TlnZ, where Z is the partition function, the 
computation for the rate boils down to the evaluation of the partition function 
for the system comprised of critical bubbles of the lower energy phase inside 
the metastable phase. 

If we imagine the space of all possible field configurations for a given 
model, there will be different paths to go from the metastable to the ground 
state. We can think of the two states as being separated by a hill of a given 
"height". The energy barrier for the decay is then related to the height 
of this hill. At the top of the hill, only one direction leads down to the 
ground state, the unstable direction. Fluctuations about this direction will 
grow, with rate given by the negative eigenvalue which appears in the above 
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formula. All other directions are positively curved, and fluctuations about 
them give rise to positive eigenvalues which do not contribute to the decay 
rate. The path which will cost less energy is the one which will dominate 
the partition function, the so-called critical bubble or bounce. It is simply 
the field configuration that interpolates between the two stable points in the 
energy landscape, the metastable and ground state. The energy barrier for 
the decay is the energy of this particular field configuration. 

For a dilute gas of bubbles only, the partition function for several bubbles 
is given by p), |§, 



Z(<p f )+Z(<p f ) 



Z(<Pb 



z (<Pf) 



+ Z(<p f )- 



Z(<p f ) 



+ 



Z(ip f )exp 



%>/)J 



(2) 



where tpf is the metastable vacuum field configuration and ^ is the bubble 
configuration, the bounce solution to the 0(3)-symmetric Euclidean equation 
of motion. We must evaluate the partition functions above. This is done 
by the saddle-point method, expanding the scalar field 0(x,r), such that 
0(x, r) -> ip f + C(x, r) for Z((p f ), and 0(x, r) -> <^ 6 (x) + 77 (x, r) for Z(<p b ), 
where £(x, r) and r/(x, r) are small fluctuations about equilibrium. 

It is crucial to note that the saddle-point, or Gaussian, method only 
gives good results if indeed the fluctuations about equilibrium are sufficiently 
small that nonlinear terms in the fields can be neglected. Even though the 
method sums over all amplitude fluctuations, it does so by assuming that 
the functional integral is well approximated by truncating the expansion 
of the action to second order. The efficiency of the method relies on the 
fact that higher amplitudes will be suppressed fast enough that their con- 
tribution to the partition function will be negligible. One can visualize this 
by comparing a sharp parabolic curve with a flatter one with minimum at 
Xq, and investigating when / dxe~^ x ^ will be well approximated by writing 
f(x) ~ f(xo) + \{x — x ) 2 f"(x ). For a sharp curve, larger amplitude fluctu- 
ations will be strongly suppressed and thus give a negligible contribution to 
the integral over all amplitudes. Clearly, this will not be the case for flatter 
curves. 

using the saddle-point method one obtains for the 
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ratio of partition functions, zfof) > 



Zyipb) saddle— point 



Z(ip f ) 



det(-n E + V^Mh 
_det(-n E + V»{<p f ))f i _ 



-AS 



(3) 



where [det(M)^]~5 = J D V exp [- j/^ dr J d 3 x^7][M}r]} and AS = S E (ip b ) - 
S E (<ff) is the difference between the Euclidean actions for the field config- 
urations ifb and iff. [Note that S E ((p), and hence AS, does not include 
any temperature corrections. It would if one had summed over other fields 
coupled to ip.] Thus, the free energy of the system is, 



-T 



det(-n E + V"(y b ))p 
det(-a E + V"(ip f )) p _ 



-AS 



(4) 



Let me stress again the assumptions that go into computing the free en- 
ergy. First, that the partition function is given by Eq. |2| within the dilute gas 
approximation, and second, that the partition function is evaluated approx- 
imately by assuming small fluctuations about the homogeneous metastable 
state iff. It is clear that for situations in which there are large amplitude 
fluctuations about the metastable equilibrium state the above formula must 
break down. Thus the breakdown of the expression for the rate is intimately 
connected with the question of how well-localized the system is about the 
metastable state as the temperature drops below the critical temperature 
T c . Homogeneous nucleation, as its name already states, is only accurate 
when the metastable state is sufficiently homogeneous. In the presence of 
inhomogeneities, there is no reason to expect that the decay rate formula 
will apply. The question then is to quantify when does it break down and 
how can we incorporate nonperturbative corrections to the decay rate in a 
consistent way. 



2 Nonperturbative corrections to decay rates 



In order to investigate the importance of large-amplitude fluctuations in the 
description of first order phase transitions, I have developed numerical sim- 
ulations in two |13| and, with J. Borrill, three ]14J spatial dimensions, which 
measured the fraction of the volume of the system in the initial phase as a 
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function of the barrier height. Since these have been documented elsewhere, 
here I quickly describe the main idea and results. 

Imagine a scalar field with a degenerate double-well potential. The field 
is coupled to a thermal bath at temperature T through a Langevin-like equa- 
tion which assumes that the bath is Markovian, i.e., the noise is white and 
additive. The system is artificially divided into two regions, left and right 
of the maximum of the potential, call it the negative and positive regions, 
respectively. The system is prepared initially in one of the regions, say, 
the negative region with 0(x, t = 0) = —1. The coupling to the bath will 
then drive fluctuations of the field around this minimum and we measure 
the fraction of the total volume in each of the two regions as a function of 
the parameters controlling the height of the potential barrier, usually the 
temperature and/or a coupling constant. 

We observed that while for large enough potential barriers the system 
remained well-localized around its initial state, a sharp change of behavior 
occurred for a critical value of the specific control parameter being varied. 
In the case examined in Ref. the control parameter was the quartic 



coupling of the scalar field, A. We showed that for A > A c the system became 
completely mixed, in that the volume was equally shared by the positive and 
negative regions. In other words, for A > A c , the effective potential describing 
the system is not a degenerate double-well, but a parabolic curve centered at 
((f)) = 0; Thermal fluctuations have "restored" the symmetry of the system. 

The challenge was thus to model the large amplitude fluctuations which 
were responsible for this phase mixing. In what follows I review the so-called 
subcritical bubbles method which can model quantitatively the dynamics of 
large, nonperturbative, thermal fluctuations in scalar field theories. 
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LECTURE II 



2.1 Modeling nonperturbative fluctuations: Symme- 
try restoration and phase mixing 

As was stressed before, the computation of decay rates based on homogeneous 
nucleation theory assumes a smooth metastable background over which crit- 
ical bubbles of the lower free energy phase will appear, grow and coalesce, 
as the phase transition evolves. However, as the results from the numerical 
simulations indicate, the assumption of smoothness is not always valid. To 
the skeptical reader, I point out that several condensed matter experiments 
indicate that homogeneous nucleation fails to describe the transition when 
the nucleation barrier (AS/T) becomes too small. Furthermore, the agree- 
ment between theory and experiment has a long and problematic history 
Homogeneous nucleation has to be used with care, in a case by case basis. 

The basic idea is that in a hot system, not only small but also large ampli- 
tude fluctuations away from equilibrium will, in principle, be present. Small 
amplitude fluctuations are perturbatively incorporated in the evaluation of 
the finite temperature effective potential, following well-known procedures. 
Large amplitude fluctuations probing the nonlinearities of the theory are not. 
Whenever they are important, the perturbative effective potential becomes 
unreliable. In an ideal world, we should be able to sum over all amplitude 
fluctuations to obtain the exact partition function of the model, and thus 
compute the thermodynamic quantities of interest. However, we can only to 
this perturbatively, and will always miss information coming from the fluc- 
tuations not included in its evaluation. If large amplitude fluctuations are 
strongly suppressed, they will not contribute to the partition function, and 
we are in good shape. But what if they are important, as argued above? 
We can try to approach this question avoiding complicated issues related 
to the evaluation of path integrals beyond the Gaussian approximation by 
obtaining a kinetic equation which describes the fraction of volume popu- 
lated by these large amplitude fluctuations. In order to keep the treatment 
simple, and thus easy to apply, several assumptions are made along the way, 
which I believe are quite sensible. In any case, the strength of the method is 
demonstrated when the results are compared with the numerical experiments 
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described before. 

Large amplitude fluctuations away from equilibrium are modelled by 
Gaussian-profile spherically-symmetric field configurations of a given size 
and amplitude. They can be thought of as being coreless bubbles. Keep- 
ing with the notation of the numerical experiment, fluctuations away from 
the 0-phase [called the "negative phase" above], and into the 0-phase are 
written respectively as, 



2 /R 2 



4o(r) = C (l - e~ r2 l R2 ) , (5) 



where R is the radial size of the configuration, and <p c is the value of the am- 
plitude at the bubble's core, away from the 0-phase. In previous treatments 
(cf. Refs. |LB| and fll7|), it was assumed that (p c = <f) + , that is, that the 



configuration interpolated between the two minima of the effective poten- 
tial, and that R = £(T), where £(T) = m(T) _1 is the mean-field correlation 
length. But in general, one should sum over all radii and amplitudes above 
given values which depend on the particular model under study. This will 
become clear as we go along. 

Define dn(R, <p, t) as the number density of bubbles of radius between 
R and R + dR at time t, with amplitudes > <p c between and <fi + d(j). 
By choosing to sum over bubbles of amplitudes <p c and larger, we are effec- 
tively describing the system as a "two-phase" system. For example, in the 
numerical simulation above it was assumed that the negative-phase was for 
amplitudes 4> < ma x, and the positive-phase was for amplitudes 4> > max - 
Clearly, for a continuous system this division is artificial. However, since the 
models we are interested in have two local minima of the free energy, this 
division becomes better justified. Fluctuations with small enough amplitude 
about the minima are already summed over in the computation of the effec- 
tive potential. It is the large amplitude ones which are of relevance here. To 
simplify the notation, from now on I will denote by "+ phase" all fluctuations 
with amplitudes <fi > C and larger. The choice of C is model-dependent, as 
will be clear when we apply this formalism to specific examples. 

The fact that the bubbles shrink will be incorporated in the time depen- 
dence for the radius Rf\. Here, I will only describe a somewhat simplified 

4 Of course, the amplitude 4> will also be time-dependent. However, its time-dependence 
is coupled to that of the radius, as recent studies have shown [|l8). In order to describe 
the effect of shrinking on the population of bubbles it is sufficient to include only the time 
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approach to the dynamics. More details are provided in the work by Gleiser, 



Heckler, and Kolb |19| . The results, however, are essentially identical. 

The net rate at which bubbles of a given radius and amplitude are created 
and destroyed is given by the kinetic equation, 

dn{R,t) dn(R,t) 
dt = 




(6) 

Here, T ^ + (R) (T + ^ (R)) is the rate per unit volume for the thermal nucle- 
ation of a bubble of radius R of positive-phase within the 0-phase (0-phase 
within the positive-phase). Vo(+) is the volume available for nucleating bub- 
bles of the +(0) phase. Thus we can write, for the total volume of the system, 
Vt = Vo + V+, expressing the fact that the system has been "divided" into 
two available phases, related to the local minima of the free energy density. 
It is convenient to define the fraction of volume in the + phase, 7, as 

I — I- 

In order to compute 7 we must sum over all bubbles of different sizes, 
shapes, and amplitudes within the + phase, i.e., starting with </> min > <f) c . 
Clearly, we cannot compute 7 exactly. But it turns out that a very good ap- 
proximation is obtained by assuming that the bubbles are spherically sym- 
metric, and with radii above a given minimum radius, -R m i n - The reason 
we claim that the approximation is good comes from comparing the results 
of this analytical approach with numerical simulations. The approximation 
starts to break down as the background becomes more and more mixed, and 
the morphology of the "bubbles" becomes increasingly more important, as 
well as other terms in the kinetic equation which were ignored. For exam- 
ple, there should be a term which accounts for bubble coalescence, which 
increases the value of 7. This term becomes important when the density of 
bubbles is high enough for the probability of two or more of them coalescing 
to be non- negligible. As we will see, by this point the mixing is already so 
pronounced that we are justified in neglecting this additional complication 



dependence of the radius. 
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to the kinetic equation. As a bonus, we will be able to solve it analytically. 
The expression for 7 is, 



roc roo / \ Q2 n 

1- / 7^ d( t )dR ■ 8 

•/tfmln J Rmin \ 3 / 9<j)dR 

The attentive reader must have by now noticed that we have a coupled 
system of equations; 7, which appears in the rate equation for the number 
density n, depends on n itself. And, to make things even worse, they both 
depend on time. Approximations are in order, if we want to proceed any 
further along an analytical approach. The first thing to do is to look for the 
equilibrium solutions, obtained by setting dn/dt = in the kinetic equation. 
In equilibrium, 7 will also be constant^. If wished, after finding the equi- 
librium solutions one can find the time-dependent solutions, as was done in 



Ref. [fT7 |. Here, we are only interested in the final equilibrium distribution 
of subcritical bubbles, as opposed to the approach to equilibrium. 

The first approximation is to take the shrinking velocity of the bubbles to 
be constant, dR/dt = —v. This is in general not the case (c/. Ref. [fL8|), but 
it does encompass the fact that subcritical bubbles shrink into oblivion. The 
strength of the thermodynamic approach is that details of how the bubbles 
disappear are unimportant, only the time-scale playing a role. The second 
approximation is to assume that the rates for creation and destruction of 
subcritical fluctuations are Boltzmann suppressed, so that we can write them 
as T = AT A e~ Fsc/T , where A is an arbitrary constant of order unity, and 
F SC (R, cf) c ) is the cost in free energy to produce a configuration of given radius 
R and core amplitude <p c . For the Gaussian ansatz we are using, F sc assumes 
the general form, F sc = aR + f3R 3 , where a = b<p 2 c (b is a combination of 
7r's and other numerical factors) and (3 depends on the particular potential 
used. In practice, the cubic term can usually be neglected, as the free energy 
of small (R ~ £) subcritical bubbles is dominated by the gradient (linear) 
term. We chose to look at the system at the critical temperature T c . For 

this temperature, the creation and destruction rates, r ^+ and r H >0 are 

identical. Also, for T c , the approximation of neglecting the cubic term is 



J This doesn't mean that thermal activity in or between the two phases is frozen; equi- 
librium is a statement of the average distribution of thermodynamical quantities. Locally 
bubbles will be created and destroyed, but always in such a way that the average value of 
n and 7 are constant. 
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very good (in fact it is better and better the larger the bubble is) even for 
large bubbles, since for degenerate vacua there is no gain (or loss) of volume 
energy for large bubbles. Finally, we use that V + /Vt = 7 in the T + ^ term. 
A more sophisticated approach is presented in Ref. |L9| . 
We can then write the equilibrium rate equation as, 



where, 



c = (1 - 2 1 )AT i /v, f(R) = e~ Fsc/T . (10) 

Integrating from -R m i n and imposing that n(R — > 00) = 0, the solution is 
easily found to be, 

n(R) = —jL- e -«*-W . (11) 
a((p c )/T 

Not surprisingly, the equilibrium number density of bubbles is Boltzmann 
suppressed. But we now must go back to 7, which is buried in the definition 
of c. We can solve for 7 perturbatively, by plugging the solution for n back 
into Eq. [8[ After a couple of fairly nasty integrals, we obtain, 

g(a((f) 

min ) 7 -^min ) ) /i o\ 

7= l + 2( 7 («(0 min ), J R min ) ' 1 ' 

where, 



4tt {AT 4 \ /T\ 3 e ~ aR ^ T 



6+(^^) 3 + 3^^(2- 



3 y v J \aj a/T 

(13) 

We can now apply this formalism to any model we wish. The first obvious 
application is to compare 7 obtained from the numerical experiment with 
the value obtained from the kinetic approach. From the definition of the 
equilibrium fractional population difference, AFeq (0 c ) = f^ q — _/+ q , 

AF EQ (fl c ) = 1 - 2 7 . (14) 
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Thus, it is straightforward to extract the value of 7 from the numerical 
simulations as a function of A. Also, as we neglected the volume contribution 
to the free energy of subcritical bubbles, we have, 

F sc = a(0 c )R min = ^n^ 2 X 2 ^e c )R min , (15) 

o 

where, as you recall, X_ is the position of the maximum of the mean-field 
potential used in the simulations. So, we must sum over all amplitudes with 
X > X_, and all radii with R > 1 (in dimensionless units), as we took the 
lattice spacing to be £ = 1. That is, we sum over all possible sizes, down to 
the minimum cut-off size of the lattice used in the simulations. In practice, 
we simply substitute <p c = X_ and -R m i n = 1 in the expression for 7. In 
Fig. 1, we compare the numerical results for 7 (dots) with the results from 
the analytical integration of the kinetic equation. The plots are for different 
values of the parameter A/ v. Up to the critical value for A ~ 0.025, the 
agreement is very convincing. As we increase A into the mixed phase region 
of the diagram, the kinetic approach underestimates the amount of volume 
in the +- phase. This is not surprising, since for these values of A the density 
of subcritical bubbles is high enough that terms not included in the equation 
become important, as I mentioned before. However, the lack of agreement 
for higher values of A is irrelevant, if we are interested in having a measure 
of the smoothness of the background; clearly, the rise in 7 is sharp enough 
that homogeneous nucleation should not be trusted for A > 0.024 or so, as 
the fraction of volume occupied by the +-phase is already around 30% of the 
total volume. Subcritical bubbles give a simple and quantitatively accurate 
picture of the degree of inhomogeneity of the background, offering a guideline 
as to when homogeneous nucleation theory can be applied with confidence, 
or, alternatively, when the effective potential needs higher order corrections. 



2.2 Modeling nonperturbative fluctuations: "Inhomo- 
geneous" nucleation 

Now we apply the subcritical bubbles method to the decay of metastable 
states in the case that the homogeneous nucleation formalism (section I) 



does not apply. Details of this work can be found in Ref. [^1[ . 

As mentioned before, if there is significant phase mixing in the back- 
ground metastable state, its free-energy density is no longer V(<ft = 0), where 



12 




0.01 0.02 0.03 0.04 0.05 
X 



Figure 1: The fraction of the volume in the + phase. The dots are from 
numerical simulations of Ref. 14, while the lines are the solutions of the 
Boltzmann eq. for different values of the parameter 



I assume the potential has a metastable state at = 0. One must also 
account for the free-energy density of the nonperturbative, large-amplitude 
fluctuations. Since there is no formal way of deriving this contribution out- 
side improved perturbative schemes, we propose to estimate the corrections 
to the background free-energy density by following another route. We start 
by writing the free energy density of the metastable state as V((f> = 0) +J- SC , 
where J- sc is the nonperturbative contribution to the free-energy density due 
to the large amplitude fluctuations, which we assume can be modelled by 
subcritical bubbles. We will calculate !F SC further below. 

We thus define the effective free-energy difference between the two min- 
ima, AV cg , which includes corrections due to phase mixing, as 

AV cg = AVo + Fsc (16) 

which is the sum of the free-energy difference calculated in the standard 
way, and the "extra" free-energy density due to the presence of subcritical 
bubbles. Henceforth, the subscript 'eg' will stand for "coarse-grained". 
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Since for degenerate potentials (temperature-dependent or not) no criti- 
cal bubbles should be nucleated, taking into account subcritical bubbles must 
lead to a change in the coarse-grained free-energy density (or potential) de- 
scribing the transition. Thus, it should be possible to translate the "extra" 
free energy available in the system due to the presence of subcritical bubbles 
in the background into a corrected potential for the scalar order parameter. 
We will write this corrected potential as V^ g (</>). 

The standard coarse-grained free energy is calculated by integrating out 
the short-wavelength modes (usually up to the correlation length) from the 
partition function of the system, and is approximated by the familiar form 

m 

F cg (<p) = J d 3 r Q(V0) 2 + K g (0)) . (17) 

How do we estimate V cg ? One way is to simply constrain it to be consistent 
with the thin wall limit. That is, as V cg (4>) approaches degeneracy (i.e. 
AV r cg (0) — > 0), it must obey the thin wall limit of eq. (|1^). Note that with 
a simple rescaling, a general polynomial potential (to fourth order) can be 
written in terms of one free parameter. Thus, the thin wall constraint can 
be used to express the corrected value of this parameter in terms of JF SC in 
appropriate units. The free energy of the critical bubble is then obtained by 
finding the bounce solution to the equation of motion V 2 — dV cg ((p) / d<f) = 
by the usual shooting method, and substituting this solution into eq. (|TT|). 

In order to determine V cg , we must first calculate the free-energy density 
J- sc of the subcritical bubbles. From the formalism presented in the previous 
subsection, 

roo /-.Rmax d 2 n( R t) 

•></>mi„ JRmin ORd(pA 

where min defines the lowest amplitude within the +phase, typically (but 
not necessarilly) taken to be the maximum of the double-well potential. -R m i n 
is the smallest radius for the subcritical bubbles, compatible with the coarse- 
graining scale. For example, it can be a lattice cut-off in numerical simula- 
tions, or the mean-field correlation length in continuum models. As for R max , 
it is natural to choose it to be the critical bubble radius. 

As an application of the above method, we investigated nucleation rates 
in the context of a 2-d model for which accurate numerical results are avail- 



able [22]. This allowed us to compare the results obtained by incorporating 
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Figure 2: Comparison between numerical data and theoretical predictions 
for the decay barrier with and without the inclusion of subcritical bubbles. 
The parameter a is related to an extra term in the Boltzmann eq. which can 
be safely neglected. 

subcritical bubbles into the calculation of the decay barrier with the results 
from the numerical simulations. The potential used was written in terms of 
one dimensionless parameter A = m 2 /i/g 2 , 

V{<1>) = \<P 2 - 1<P 3 + ^\ (19) 

This double-well potential is degenerate when A = 1/3, and the second min- 
imum is lower than the first when A < 1/3. 

As argued before, we find the new coarse-grained potential V cg (or, equiv- 
alently, A cg ) by constraining it to agree with the thin wall limit. Simple 
algebra from eqs. flTBD and (|T9|) yields, to first order in the deviation from 
degeneracy, 

A cg = A - ^ (20) 

2 

where JF SC = - 2 -g-^ r sc is the dimensionless free-energy density in subcritical 
bubbles. The new potential V cg is then used to find the bounce solution and 
the free energy of the critical bubble. 
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In Fig. 2 we show that the calculation of the nucleation barrier including 
the effects of subcritical bubbles is consistent with data from lattice simula- 
tions, whereas the standard calculation overestimates the barrier by a large 
margin. In fact, the inclusion of subcritical bubbles provides a reasonable 
explanation for the anomalously high nucleation rates observed in the simu- 
lations close to degeneracy. 



3 Matching numerical simulations to contin- 
uum field theories 

As a final topic to be discussed in this lecture, I would like to change gears 
and briefly turn to the issue of how to match numerical simulations of field 
theories with their continuum counterparts. In particular, I am interested 
in situations where a degree of stochasticity is present in the simulations, 
as for example happens when we model the coupling of fields to a thermal 
or quantum bath via a Langevin-like equation, or even in the form of noisy 
initial conditions. 

Although field theories are continuous and usually formulated in an in- 
finite volume, lattice simulations are discrete and finite, imposing both a 
maximum ( "size of the box" L) and a minimum (lattice spacing 8x) wave- 
length that can be probed by the simulation. When the system is coupled 
to an external thermal (or quantum) bath, fluctuations will be constrained 
within the allowed window of wavelengths, leading to discrepancies between 
the continuum formulation of the theory and its lattice simulations; the re- 
sults will be dependent on the choice of lattice spacing. 

Parisi suggested that if proper counterterms were used, this depedence on 
lattice spacing could be attenuated p3| . Recently, Borrill and Gleiser (BG) 
have examined this question within the context of 2-d critical phenomena 



2J]. They have computed the counterterms needed to render the simu- 
lations indepedent of lattice spacing and have obtained a match between 
the simulations and the continuum field theory, valid within the one-loop 
approximation used in their approach. Here, I want to focus mostly on the 
application of these techniques to 1-d field theories, in particular to the prob- 
lem of thermal nucleation of kink-antikink pairs. [This is based on work of 



Ref. p|] 
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Even though 1-d field theories are free of ultra-violet divergences, the 
ultra-violet cutoff imposed by the lattice spacing will generate a finite contri- 
bution to the effective potential which must be taken into account if we are 
to obtain a proper match between the theory and its numerical simulation 
on a discrete lattice. If neglected, this contribution may compromise the 
measurement of physical quantities such as the density of kink-antikink pairs 
or the effective kink mass. 

For classical, 1-dimensional finite-temperature field theories, the one-loop 



corrected effective potential is given by the momentum integral |23 



T r°° dk 



1 + 



k 2 



+ ( 21 ) 



As mentioned before, the lattice spacing 5x and the lattice size L in- 
troduce long and short momentum cutoffs A = n/5x and k min = 2n/L, 
respectively. Lattice simulations are characterized by one dimensionless pa- 
rameter, the number of degrees of freedom N = L/Sx. For sufficiently large 
L one can neglect the effect of k min and integrate from to A. For V " <C A 2 
(satisfied for sufficiently large A), the result can be expanded into 



Kl(0,A) = V + 



T i — T V" 

-\Ao" r~ + AT O 

4 v 4vr A 



A 4 



(22) 



As is to be expected for a 1-dimensional system, the limit A — > oo exists 
and is well-behaved; there is no need for renormalization of ultra-violet diver- 
gences. However, the effective one-loop potential is lattice-spacing dependent 
through the explicit appearance of A, and so are the corresponding numerical 
simulations. In order to remove this dependence on 5x, we follow the renor- 
malization procedure given by BG p3}| ; it is irrelevant if the A-dependent 
terms are ultra-violet finite (d = 1) or infinite (d > 2). In the lattice formu- 
lation of the theory, we add a (finite) counterterm to the tree-level potential 
Vq to remove the lattice-spacing dependence of the results, 



t TO) 

4:71 A 



(23) 



There is an additional, A-independent, counterterm which was set to zero 
by an appropriate choice of renormalization scale. The lattice simulation then 
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t t 

Figure 3: Average field value <fi(t) for T = 0.1 using the tree-level potential, 
left, and the corrected potential, right. The filter cutoff is AL = 3. 



uses the corrected potential 

TSx 

Viatt(0) = V (<f>) + ^T y o'(0) • (24) 

Note that the above treatment yields two novel results. First, that the 
use of Vi att instead of V gets rid of the dependence of simulations on lattice 
spacing. [Of course, as 5x — > 0, V^tt "~ > V . However, this limit is often not 
computationally efficient.] Second, that the effective interactions that are 
simulated must be compared to the one-loop corrected potential Vil(0) of 
Eq. (pl|); once the lattice formulation is made independent of lattice spacing 
by the addition of the proper counterterm(s), it simulates, within its domain 
of validity, the thermally corrected one- loop effective potential. 

Applying this method to the formation of kink-antikink pairs, we get a 
corrected potential, 

Matt(0) = Vo((/>) + ^ AT( ^ 2 ! (25) 
simulations using V[ at t will, in principle, match the continuum theory 

V 1L te) = Vte) + ^VV-^O , (26) 
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which has (shifted) minima at ±0 m i n (T), with m i n (T) < 0o- 

As a first test of our procedure, we investigate the mean field value (j)(t) = 
(1/L) / <p(x,t)dx before the nucleation of a kink-antikink pair, i.e., while the 
field is still well localized in the bottom of the well. In Fig. 3 we show 
the ensemble average of (after 100 experiments) for different values of Sx, 
ranging from 1 down to 0.1, at T = 0.1. The simulations leading to the left 
graphs use the "bare" potential Vo, whereas the right graphs are produced 
employing Vi att (Eq. ||). 

Perhaps the most difficult task when counting the number of kink-antikink 
pairs that emerge during a simulation is the identification of what precisely 
is a kink-antikink pair at different temperatures. Typically, we can identify 
three "types" of fluctuations: i) small amplitude, perturbative fluctuations 
about one of the two minima of the potential; ii) full-blown kink-antikink 
pairs interpolating between the two minima of the potential; iii) nonpertur- 
bative fluctuations which have large amplitude but not quite large enough 
to satisfy the boundary conditions required for a kink-antikink pair. These 
latter fluctuations are usually dealt with by a smearing of the field over a cer- 
tain length scale. Basically, one chooses a given smearing length AL which 
will be large enough to "iron out" these "undesirable" fluctuations but not 
too large that actual kink-antikink pairs are also ironed-out. The choice of 
AL is, in a sense, more an art than a science, given our ignorance of how to 
handle these nonperturbative fluctuations. 

The smearing was implemented as a low pass filter with filtering cut- 
off AL; the field is Fourier transformed, filtered at a given wavelength, and 
Fourier transformed back. We counted pairs by identifying the zeros of the 
filtered field. Choosing the filter cutoff length to be too large may actually 
undercount the number of pairs. Choosing it too low may include nonper- 
tubative fluctuations as pairs. We chose AL = 3 in the present work, as 
this is the smallest "size" for a kink-antikink pair. In contrast, in the works 
by Alexander et al. a different method was adopted, that looked for zero- 
crossings for eight lattice units (they used 5x = 0.5) to the left and right 
of a zero crossing We have checked that our simulations reproduce 



the results of Alexander et al. if we: i) use the bare potential in the lat- 
tice simulations and ii) use a large filter cutoff length AL. Specifically, the 
number of pairs found with the bare potential for T = 0.2, Sx = 0.5 are: 
n p = 36, 30, and 27, for AL = 3, 5, and 7 respectively. Alexander et al. 
found (for our lattice length) n p = 25. Comparing results for different AL, 
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Figure 4: Density of kink-antikinks (half of density of zeros), for T = 0.2 and 
5x = 1, 0.5, 0.2, and 0.1. The filter cutoff is AL = 3. 

it is clear that the differences between our results and those of Alexander et 
al. come from using a different potential in the simulations, viz. a dressed 
vs. a bare potential. For small 5x these differences disappear. 

Fig. 4 compares measurements of the kink-antikink pair density (half the 
number of zeros of the smeared field), ensemble-averaged over 100 experi- 
ments, for different lattice spacings. Again it is clear from the graphs on 
the left that using the tree-level potential Vq in the simulations causes the 
results to be dependent on 5x, whereas the addition of the finite counterterm 
removes this problem quite efficiently 

Another step is to establish what is the continuum theory being simulated. 
Due to space limitations, I refer the reader to the work of Ref. PS| for more 
details. 

This work was written in part while the author was visiting the Osser- 
vatorio Astronomico di Roma and the Nasa/Fermilab Astrophysics Center. 
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